Magneto-transmission of multi-layer epitaxial graphene and bulk graphite: 

A comparison 
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Magneto-transmission of a thin layer of bulk graphite is compared with spectra taken on multi- 
layer epitaxial graphene prepared by thermal decomposition of a SiC crystal. We focus on the 
spectral features evolving as vB, which are evidence for the presence of Dirac fermions in both 
materials. Whereas the results on multi-layer epitaxial graphene can be interpreted within the 
model of 2D Dirac fermions, the data obtained on bulk graphite can only be explained taking into 
account the 3D nature of graphite, e.g. by using the standard Slonczewski-Weiss-McClure model. 

PACS numbers: 71.70.Di, 76.40.+b, 78.30.-j, 81.05.Uw 



I. INTRODUCTION 

Within the last few years, the properties of massless 
and massive "relativistic" particles have been probed 
in magneto-optical experiments carried out on sev- 
eral carbon-based systems: graphene monolayer^ and 
bilayer, 3 multi-layer epitaxial graphene as well as 
bulk graphite i 8 ' 9 ' 10 ! 11 These experiments offer an impor- 
tant insight into the electronic structure of these materi- 
als, in particular, they directly demonstrate the unusual 
\/"B-scaling of Landau levels (LLs) of Dirac fermions. 
Naturally, these experiments stimulated intensive theo- 
retical work in this areai 12 ' 13 i 14 i 15 ' 16 

Recently, remarkable experiments have been per- 
formed on the same systems without applied magnetic 
field. A universal value of optical conductivity was 
demonstrated in graphene 1 ^ as well as in bulk graphite. 18 
The optical response of graphene monolayer and bilayer 
samples as a function of the gate voltage was reported 
in several experimental works, see Refs. [l3l2Gl2lll22l . and 
this topic was addressed also theoreticallyi 23 ' 24 ^ First 
results of time-resolved optical experiments on multi- 
layer graphene are also availablej 26 ' 27 

In this paper, we present a direct comparison of the 
magneto-transmission spectra taken on a thin layer of 
bulk graphite and multi-layer epitaxial graphene. We 
focus on spectral features exhibiting a v^B-dependence, 
which corresponds to the optical response of Dirac 
fermions, and show that there are essential differences 
in the spectra of both materials. 



II. EXPERIMENTAL DETAILS 

The investigated multi-layer graphene sample was pre- 
pared by thermal decomposition on the carbon face 
of a 4H-SiC substrat o 28 i 29 and contains around ^100 



graphene layers. The thin layer of bulk graphite was pre- 
pared by simple exfoliation as described in Ref. As no 
significant difference in FIR optical response of natural 
graphite and HOPG was found^ we present only results 
taken on natural graphite due to its higher crystalline 
quality. Both samples were characterized using micro- 
Raman. In the case of multi-layer epitaxial graphene, 
the Raman signature of decoupled layers, equivalent to a 
single flake of exfoliated graphene, 30 as well as of addi- 
tional graphite residuals were found, see Ref. The 
Raman spectra taken on a thin graphite layer domi- 
nantly showed a multi-component 2D band typical of 
many Bernal-stacked sheets, but some minor traces of 
decoupled layers were found, see the discussion in Ref. [ToL 
To measure the FIR transmittance of the sample, the 
radiation of globar, delivered via light-pipe optics to the 
sample and detected by a Si bolometer placed directly 
below the sample, was analyzed by a Fourier transform 
spectrometer. All measurements were performed in the 
Faraday configuration with the magnetic field applied 
normal to graphenc/graphitc layers. All the spectra were 
taken with non-polarized light in the spectral range of 5- 
350 meV, limited further by several regions of low tape 
transmissivity or the SiC opacity, see grey areas in Fig.[T] 



III. RESULTS AND DISCUSSION 

A comparison of the transmission spectra taken on 
multi-layer epitaxial graphene and a thin layer of bulk 
graphite is presented in Fig. [TJ Starting with results 
taken on multi-layer layer graphene in Fig. [T£i, a series 
of absorption lines is observed and denoted by Roman 
letters, following the notation introduced in Ref. 1. As- 
suming the LL spectrum of graphene, E n — cy2ehBn, 
together with the corresponding selection rules for dipolc- 
allowed transitions An = ±1, the observed absorption 
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lines B,C,D,E and F can be clearly identified as inter- 
LL transitions L_ m — *L m _|_i and L_( m+ x) ^L m with 
m = 0,1,2,3 and 4, respectively. The Fermi velocity 
is found to be c = (1.02 ± 0.01) x 10 6 m.s^ 1 . Hence, 
these results are fully consistent with a model of 2D Dirac 
fcrmions and individual sheets in multi-layer epitaxial 
graphene indeed behave as if they are electronically de- 
coupled. This was recently explained by a random mu- 
tual rotation of adjacent sheets in this material* 3 ^ Similar 
spectra were also measured on a single sheet of exfoli- 
ated graphene, 1 however, the Fermi velocity seems to be 
enhanced by amount of «10% and these data also sug- 
gest some influence of electron-electron interaction, not 
observed in transmission spectra of multi-layer epitaxial 
graphene. 

The results obtained on a thin layer of bulk graphite, 
see Fig. lb, are more complex. Basically, all the ab- 
sorption lines observed in multi-layer graphene arc also 
found in spectra of bulk graphite with practically the 
same Fermi velocity c = (1.02 ± 0.02) x 10 6 m.s -1 . This 
justifies the same notation using the Roman letters. In 
addition to these lines, another series of transitions, de- 
noted by Greek letters and clearly exhibiting the \f~B- 
dependence, is present in transmission spectra. These 
lines cannot be assigned to any dipole-allowed transitions 
between LLs in graphene. Nevertheless, the energies of 
the additional lines a, 7, S and e exactly match to tran- 
sitions symmetric around the Dirac point, L_ m — >L m , 
with indices m = 1,2,3 and 4, respectively. The (3 line 
can be identified as transitions L_!(_ 3 ) — >L 3 ( 1 ). 

To explain the absorption lines denoted by the Greek 
letters we have to abandon the simplified model of 
2D Dirac fermions and consider the full band struc- 
ture of bulk graphite. According to the standard SWM 
model , 33 ' 34 Dirac fermions are located in the vicinity of 
the H point, where the bands E\, E 2 and the doubly de- 
generate E3 are close to the Fermi level. If the magnetic 
field is applied along the c-axis of the crystal, Landau 
levels or more accurately, Landau bands are created, see 
Fig. [21 having at the H point (k z = 0.5) an analytical 
foimM (n > 1): 

Kf 1 = 
El = A 

£§± = E$ 1 = ~ ± (1) 

where A denotes the pseudogap, i.e. the distance of E\ 
(and also E2) and E3 bands at k z = 0.5 and B = 0. The 
parameter £ is related to the Fermi velocity as £ = 2c 2 eh. 

The form of LLs {1} implies several important H point- 
related optical properties of graphite. It suggests that the 
LL energy spectrum typical of graphene is present also 
at the H point in graphite, when the pseudogap |A| is 
small in comparison to the energies of the LLs. Experi- 
mentally,! A[ is found to be significantly below 10 meV, 
see Refe.la3a37l and therefore, the magnetic fields above 
100 mT are sufficient to insure this condition. On the 
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FIG. 1: Transmissions spectra of multi-layer epitaxial 
graphene (a) and bulk graphite (b) for selected magnetic fields 
at T — 2 K. The absorption lines corresponding to dipole- 
allowed transitions in graphene are denoted by Roman letters. 
Greek letters are used for additional transitions which scale 
as vB and are only found in spectra taken on bulk graphite. 
For clarity, the spectra in part (a) and (b) were shifted by 
0.18 and 0.10, respectively. 

other hand, each LL in Eq. (JTJ) is (with the exception 
of E° and E^ 1 ) doubly degenerate, Eg + = E^ 1 and 
i?3_ = E$ , see Fig. [2j This double degeneracy is 
in addition to the spin and valley degeneracies present 
in graphene. Obviously, taking account of the same se- 
lection rules An = ±1, we obtain a considerably richer 
set of possible dipole-allowed transitions for graphite in 
comparison to graphene. For instance, the transition 
L_i — >Li is strictly forbidden in graphene, neverthe- 
less the absorption line at this energy (a transition) is 



3 




0.46 0.47 0.48 0.49 0.50 0.47 0.48 0.49 0.50 0.47 0.48 0.49 0.50 



(a) Momentum k z (rel. u.) (b) Momentum k z (rel. u.) ( c ) Momentum k z (rel. u.) 

FIG. 2: LL structure of bulk graphite in the vicinity of the H point. Dipole- allowed transitions between LLs depicted in blue, 
dark yellow and red colors in part (a), (b) and (c) contribute to the B,q and C lines in the transmission spectrum, respectively, 
cf. Fig. [Tp. The parameters A = —5 meV and c = 1.02 x 10 6 m.s -1 are used for this calculation. The effects of trigonal 
warping are neglected. 



observed in bulk graphite at the H point due to the 
dipole-allowed transitions E\_ — > E® and E% — > -E31, 
see Fig. [2Jd. Similarly, the lines 7,(5 and e are detected 
due to transitions E^+ l -> E™ and Ef -> Eg* 1 for 
m = 1,2 and 3, respectively. The remaining f3 line can 
be identified as E\_ — > E\ and E\ — > E\+. Hence, the 
presence of absorption lines denoted by Greek letters is 
qualitatively consistent with the standard SWM model of 
graphite. Nevertheless, we are aware of inconsistency of 
our results with predictions of Koshino and Ando.— They 
expect that for graphite layer with several tens or more 
Bernal-stacked sheets, the intensities of V^-dependent 
absorption lines should be negligible in comparison with 
spectral features arising in massive electrons around the 
K point, which evolve nearly linearly with B. 

Another difference of spectra Figs. QJi and b is an ap- 
parent asymmetry of the B line in bulk graphite at low 
magnetic fields, whereas the equivalent line in the spec- 
tra of multi-layer graphene remains perfectly symmetric. 
To explain this difference, we must extend our consider- 
ations to the near vicinity of the H point and take into 
account the fcz- dependence of LLs. Such LL structure 
calculated within the standard SWM model is shown in 
Fig. UK where the two components, Bo and B_i of which 
the B line consists, are denoted by arrows. Note that the 
relatively low value of the pseudogap |A| does not allow 
to directly resolve the splitting of the B line into these 
components. At low magnetic field, when the Landau 
band E\_ is nearly parallel with E® and/or when it is 
partially or completely located above the Fermi level, we 
can expect rather asymmetrical shape of the B line with 
a well-pronounced high energy tail due to the transition 



E® — ► E§ . With increasing magnetic field, this asym- 
metry should gradually disappear, as indeed observed ex- 
perimentally, see Fig.[T|3. From the lineshape, we roughly 
estimate the position of the Fermi level slightly below ss- 
20 meV, but the situation is here complicated by the fact 
that also decoupled graphene layers in graphite, revealed 
by the Raman spectroscopy or detected in scanning tun- 
neling spectroscopy experiments^ can contribute to the 
final shape of the B line. 

Note, that for the sample of multi-layer graphene in- 
vestigated here, the fully symmetrical B line survives 
in transmission spectrum down to magnetic fields of 
40 mT, when the Fermi level is located at ss7 meV 
from the Dirac points The presented transmission exper- 
iment thus probes mainly undoped graphene sheets (no ~ 
5 x 10 9 cm -2 ), which are located further from the SiC 
substrate. The layer(s) in the immediate vicinity of the 
substrate are highly conductive (above no ~ 10 12 cm -2 ) 
and can be probed in transport experiments j ' 

Analogous considerations based on the shapes of Lan- 
dau bands in the vicinity of the H point, can also qual- 
itatively explain the significant differences in widths of 
individual absorption lines. For instance, the C line is 
roughly three times broader in comparison to the a line, 
see Fig.QJi, in spite of the fact that they both originate in 
the same LLs. For a simple explanation of this fact, wc 
just need to realize that the a line is composed of tran- 
sitions E\_ -> El and E% -> E\ + . If we take account 
of the near vicinity of k z — 0, we find out that both 
pairs of these LLs are nearly parallel, see Fig. [2b, and 
give thus a rather sharp a line in comparison to the C 
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line. This line is composed of transitions — » E^£ > 

1(2") 2(1") 

and E 2 y 1 -> E{ y \ see Fig. [2}:, which occur between 
pairs of LLs having their derivatives of opposite signs, 
which causes the broadening of the C line. This expla- 
nation can be straightforwardly generalized for all tran- 
sitions, giving qualitative explanation of the differences 
in linewidths simply using the standard SWM model. 

In principle, another and even more appealing inter- 
pretation of our data could assume that we actually 
see the response of decoupled graphene layers instead 
of the H point of bulk graphite. The presence of ad- 
ditional "Greek" lines in spectra could be then a con- 
sequence of some perturbation of the graphene system 
leading to a relaxation of the selection rules An = ±1. 
Even though the appearance of some decoupled graphene 
sheets in our sample is very likely, we have the follow- 
ing arguments to support our "bulk" interpretation: i) 
Our experiments^^ were performed on several samples 
prepared from two different types of bulk graphite and 
the almost identical results, concerning the presence, po- 
sition, lineshape and mainly mutual intensities of all 
the \/i?-scalcd absorption lines were obtained, ii) The 
"Greek" transitions show significantly smaller linewidth 
in comparison to "Roman" ones, e.g. compare lines a 
and C, which is unlikely to be caused by any perturba- 
tion relaxing the selection rules in graphene. 

In both materials, spectral features evolving nearly lin- 
early with B are also observed, see low energy part of 
spectra in Figs. [1^, and b, and represent thus an evi- 
dence for massive particles. In bulk graphite, can we 



relate these transitions to the K point, i.e. to massive 
electrons, as discussed already in Ref. 0. The optical 
response of massive electrons in the magnetic field was 
thoroughly analyzed e.g. in Refs. dHH. In case of multi- 
layer graphene, the spectral features linear in B cannot 
be explained unless some part of the sample is created 
from bulk graphite or at least from few-layer graphene 
stacks. Indeed, the traces of Bernal-stacked layers were 
found in micro-Raman experiment, as mentioned above. 



IV. CONCLUSIONS 

We have compared the optical response of Dirac 
fcrmions in bulk graphite and multi-layer epitaxial 
graphene in low magnetic fields. Whereas the results ob- 
tained on multi-layer graphene fully correspond to expec- 
tations for dipole-allowed transitions in a 2D gas of Dirac 
particles, the transmission spectra taken on bulk graphite 
appear to be more complex. The standard Slonczewski- 
Weiss-McClure model is found, in the latter case to be 
sufficient to explain the existence of all absorption lines 
scaling as \f~B~ as well as their individual lineshapes. 
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